Inflation and Reheating in Spontaneously Generated Gravity 
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Inflation is studied in the context of induced gravity (IG) 70- R, where R is the Ricci scalar, a 
a scalar field and 7 a dimensionless constant, and diverse symmetry-breaking potentials V{g) are 
considered. In particular we compared the predictions for Landau-Ginzburg (LG) and Coleman- 
Weinberg (CW) type potentials and their possible generalizations with the most recent data. We find 
that large field inflation generally leads to fewer constraints on the parameters and the shape of the 
potential whereas small field infiation is more problematic and, if viable, implies more constraints, 
in particular on the parameter 7. We also examined the reheating phase and obtained an accurate 
analytical solution for the dynamics of inflaton and the Hubble parameter by using a multiple scale 
analysis (MSA). The solutions were then used to study the average expansion of the Universe, the 
average equation of state for the scalar field and both the perturbative and resonant decays of the 
inflaton field. 

PACS numbers: 98.80Cq 



I. INTRODUCTION 



Many years ago a model for a varying gravitational 
coupling was introduced [l|. The model consisted of a 
masslcss scalar field whose inverse was associated with 
the gravitational coupling. Such a field evolved dynam- 
ically in the presence of matter and led to cosmologi- 
cal predictions differing from Einstein Gravity (EG) in 
that one generally obtained a power-law time dependence 
for the gravitational coupling. Subsequently it was sug- 
gested that the gravitational constant is generated as a 
one loop effect in some fundamental interaction [2| or 
through spontaneous symmetry breaking (for reviews see 
[3, |j]). In order to reduce the strong time dependence 
in a cosmological setting, which remained after the in- 
troduction of matter, a simple globally scale invariant 
model for induced gravity (IG) involving a scalar field 
(7 and a quartic potential Act^ was introduced [5|. The 
spontaneous breaking of scale invariance in such a con- 
text, either through the presence of a condensate [a-Q or 
quantum effects [§, [i3| , then led to EG plus a cosmolog- 
ical constant and, on treating matter as a perturbation, 
a time dependence for the scalar field (gravitational con- 
staiit) and consistent results Q. A more detailed analy- 
sis [ll[ of such a simple model including both radiation 
and matter showed that it led to EG plus a cosmological 
constant as a stable attractor among homogeneous cos- 
mologies and was therefore a viable Dark Energy model 
for a range of scalar field initial conditions and a posi- 
tive 7 coupling to the Ricci scalar ^a^R. In that earlier 
study we considered values for the scalar field which were 
sufficiently close to the spontaneously broken symmetry 
equilibrium values for the scalar field and compared our 
results with present values of the cosmological constant 
and the solar system data. In a later note [12] we studied 
the above approach for values of our parameters suffi- 



ciently far from the equilibrium values (back in time) for 
sufficient inflation and the subsequent reheating to take 
place and examined the compatibility of the predictions 
for different symmetry breaking potentials with the cur- 
rent data. 

The aim of this paper is to investigate in detail the dy- 
namics of the scalar field and the generation of scalar and 
tensor perturbations in the early Universe. 
The paper is organized as follows: in section II the for- 
malism and the general equations for the homogeneous 
dynamics are described. In section III we formalize the 
problem of perturbations in the IG context and section 
IV is dedicated to comparing the predictions of diverse 
symmetry-breaking potentials with the most recent data. 
In section V we apply a Multiple Scale Analysis (MSA) 
to the phase of coherent oscillations of the inflaton field 
at the end of inflation. In section VI we study the per- 
turbative decay of the inflaton field in the IG context 
and in section VII the resonant decay of the inflaton is 
discussed. Finally, in section VIII, our conclusions are 
summarized. 



II. IG INFLATION 

We consider the system described by the action 



S ~ I d x^—g 



^i^df,ad,a + ^a'R-V{a) 



2 



(1) 



where 7 is a dimensionless, positive definite parameter 
giving the non-minimal coupling between the scalar field 
and gravity. The Einstein-Hillbert term for gravity is re- 
placed by an effective mass-like term for the scalar field 
(T. In turn such a scalar field is also responsible for the ex- 
istence of space-time curvature and its expectation value 
dynamically generates an effective Planck mass. If we 



restrict our analysis to the homogeneous dynamics and 
we assume a spatially flat Robertson- Walker background 



rfs2 = g^^dx^'dx" = -dt^ + a^{t)dx^ 



(2) 



the variation of the above Lagrangian leads to the follow- 
ing set of independent equations 
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where we defined Kff,cr = dV/da — iV/a. The l.h.s. of 
Eq. ([3]) is positive definite or zero whereas the l.h.s of the 
same equation can also be negative in some regions of the 
phase-space. As a consequence one has some restrictions 
on the allowed phase-space. In particular one needs 



67 - v/67 (1 + 67) < -— < 67 + ^67(1 + 67) (5) 

in order for Q t o be well defined. 

It is known |13l . ll4| that by a conformal transformation: 



_ q2 



da^ 



V = vi-'^v 



^' + '^^da^ 



(6) 



where ft^ = 7cr2/Mp , wc can rewrite (up to a boundary 
term) the action in Eq. ([T]) as: 



Sf 



d^x 



— a^cra^^cr -I- V{cr) 



(7) 

and we have introduced a (reduced) Planck mass Mp 
(= (87rG)"^/2 ^[i]^ Q the Newton constant). 

The variation of Eq. ([7]) then leads to a set of equations 
analogous to those obtained from Eq. ([T]). Let us note 
that the spectrum of curvature perturbations and the 
amplitude of gravitational waves obtained are invariant 
under conformal transformations [15| : it is for this reason 
that inflationary calculations are often performed in the 
Einstein frame. 

We observe that other important quantities in cosmology 
are not left invariant under conformal transformation. 
This is the case for the Hubble parameter H. If we are 
interested in late time cosmology and use observational 
data to constrain iJ in a scalar-tensor theory [Il|, this 
should be done for the Hubble parameter in the Jordan 
frame, which is different from the Hubble parameter in 
the Einstein frame. To conclude, for us, the Jordan frame 
is the physical one. 



A. Hubble and scalar field flow functions 

As emphasized in [l2|, cosmological linear perturba- 
tions depend not only on the derivatives of the Hubble 



parameter, but also on the derivative of the scalar field 
itself. It is therefore useful to introduce the hierarchy of 
scalar field flow functions (5„ {d\n\Sn\/dN = Sn+i with 
n > 0, Sq = a/a{ti)) in addition to e„ {dlii\en\/dN = 
£)i+ij £0 = H{ti)/H) with n > 0, where ti is some initial 



time and A^ 



In 



a{U} 



is the number of e-folds. These 



two hierarchies are related by: 
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Si f Si 



l + 5i V27 



261 + S2-I 



(8) 



The above parameters are related to the Hubble flow 
functions in the Einstein frame (e^) by 



ei 



£2 



(l + 67)<52 
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(9) 



These hierarchies arise naturally with the slow-roll (SR) 
of the inflaton field. Note that in EG the equations gov- 
erning the dynamics of the scalar and tensor fluctuations 
during inflation can be written in terms of the Hubble 
flow function hierarchy. In IG the equivalent set of equa- 
tions cannot be written only in terms of Hubble flow 
function hierarchy and the scalar field flow function hi- 
erarchy is also needed. 



Homogeneous Dynamics as flow of the Hubble 
and Scalar Field Flow^ Functions 



The homogeneous dynamics of the field-gravity sys- 
tem and the slow-roll conditions for inflation have a pe- 
culiar role in the theory of cosmological perturbations 
since they provide an approximate method to determine 
the dynamics of these perturbations and compare the- 
oretical models with observations. Exact solutions for 
such a dynamics can be also found for particular choices 
of the inflaton potential both in the EG and IG frame- 
work. On writing the equations of motion in terms of the 
SR parameters e^ and Si one can easily find a set of these 
solutions. In EG one has 



Si 



Mp^ 



= ci =^ di + d2 = — 



(10) 



and the Klein-Gordon equation for the scalar field can be 
rewritten as 



■Si-h 



SidlnV 
3+— — — ^(3-ei) 
ei diner 



0. 



(11) 



From Eqs. (|10llip one easily observes that no solu- 
tion with 5i and ei simultaneously constant and dif- 
ferent from zero exists while a non trivial solution can 
be found for the case £2 = 0, (5i = zb-^/eiMp/o- and 
V ex exp (-y/eicr/Mp) , namely the well-know power-law 
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FIG. 1: The figure represents ei in Ifl3\) as a function of n 
and 7. The lighter grey region is for ei > 1, the darker grey 
region is for ei < and the intermediate grey region is for 
< ei < 1. The lines for n = 2,4 represents the de Sitter 
solutions. The darkest grey areas are the domains allowed by 
the present experimental constraints \26\) . 



inflationary solution. 

In the IG context Eqs. (jlOllip are replaced by Eq. 

and 



ei 



3 Ui - 47 + 7 



+ 5i ((5i + 52) 



dlnF 
din cr 
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5i — 67 
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Despite their quite involved form one can still see that an 
exact, non trivial solution exists with 62 = (^2 = when 
y = Vo (cr/cTo)", that is 



<5i 



7 (n - 4) 



:. ei 



7(ri-2)(n-4) 



l + 7(n + 2)' 2 + 27(^ + 2) 

On solving Eqs. (|8ll2p one obtains 



(13) 



which are not compatible with phase-space constraints. 
In Fig. ([1]) we plotted the behavior of ei in P^ on 
varying 7 and n. Decelerated, accelerated and super- 
accelerated solutions can be found inside the boundaries 
of phase space as can be seen in Fig. ([ij where the lighter 
region represents decelerated trajectories (ei > 1), the 
darker area is for super-accelerated solutions (ei < 0) 
and the intermediate grey region represents the accel- 
erated trajectories (0 < ei < 1). The outer bound- 
aries of the lighter gray area are those of phase space 
and are also those of the region of stability for such ex- 
act solutions with respect to homogeneous perturbations. 
The lines n = 2 and n — A dividing accelerated and 
super-accelerated trajectories represents the De Sitter so- 
lutions. Note finally that the exact solutions with n > 4 
have (5i < (Large Field dynamics) while for n < 4 one 
has (5i > (Small Field dynamics). 



III. COSMOLOGICAL PERTURBATIONS 

Scalar curvature perturbations produced by quantum 
fluctuations of the infiaton during the accelerated stage 
are described by ^^.(a;) = —HSa{x)/& in the uniform cur- 
vature gauge [l6|, where Sa{x) is the scalar inflaton per- 
turbation and is the correct field variable to quantize. 
The Fourier component 6ak {"R-k = —H5crk/a) of the in- 
flaton fluctuation, in the IG context, has been shown to 
satisfy the following differential equation |16j : 



5ak + \ZH+-\5<Jk 



e 



a? a^ZaSi 



[a^Z {a6i 



where 



H^a^{l + Q-i) 1 + 67 



Sak^Q (15) 



(16) 



{a + HaY (l + (5i)2- 

The above equation can be rewritten for Sk = avZ Sak 



as: 



d^5fc 



+ [fc2 + mI{tj)] Sk=0 



(17) 



(5i = 67 ± ^67(1 + 67), ei = 3 + 2Si (14) where 



Mliy) = -H^ 



Si + Sl + {3- ei)(<5i +62 + 1) + S2S3 + 



5iS2 
1 + Si 



r 



ei + 61 — 3^2 — ^3 + 



2S1S2 
1 + 61 



-2-1 



(18) 



with ?7 the conformal time (a(77)d77 = dt) and "H — inflation. In IG the Fourier modes of tensor perturba- 
a^^da/dr]. Gravitational waves are also produced during 



tions satisfy the following equation: 



hs,k + {'SH + 2H6i)hs,k + ^hs,k^0 (19) 



where s = + , x denotes the two polarization states. On 
setting Tg^k = -j^acTy/^hs.k the above equation can be 
rewritten as: 






+ [k^ + A4{7])]Ts,k = 0, (20) 



where 

M^{ri) = -n^ [2 - ei + <5i(3 + Si + S2 - ei)] . (21) 

We define the power spectra of scalar curvature pertur- 
bations and tensor perturbation as 



rnik)^^\nk\'^rn{h)(^ 



(22) 



and 



9k^ / A- ^ "* 



respectively, where fc* is a suitable pivot scale. 



(23) 



A. Exact Solutions for V{a) oc a" 

In IG power-law potentials for a lead to power-law 
inflation [l2| and the spectral index for such background 
solutions is: 



2-f{n-4y 



rig — 1 = nt = - 

7 (n- 4)^ -2 (67 + 1) 

and the exact tensor-to-scalar ratio is given by 

Vh{k) 8nt 



Vnik) 



1 "t 
^ 2 



(24) 



(25) 



which agrees with the consistency condition of power-law 
inflation in EG. 

From the most recent compilation of data [l7| we ob- 
tain the constraint: 



67 -M 
39 < , ^^,„ < 123 
7(n- 4)2 



(26) 



at the 95 % confidence level which is independent of 7 
for 7 large. 



B. Slow Roll Results 

For generic potentials the form of the exact solutions 
of Eqs. (|8I12[) is unknown and one must employ the SR 
approximation in order to obtain analytical estimates for 



the spectra of perturbations. In such an approximation 
Eqs. pi4p become 
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(27) 
(28) 



On assuming Si <^ 1 (and consequently e^ <C 1) one can 
thus neglect their variation in time and rely, to first order 
in the SR parameters, on the results obtained for exact 
evolution with S2 ^ Si y^ 0. Such an approach leads to 
the following expressions for the spectra: 



where 
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4^2(1 + 67)^2^^2 
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C = 2(2 
and 



In 2 — 5), fe is the Euler-Mascheroni constant 
2iA^CS2.)H^ 



Vhik,) 



7^2^0-2 



The spectral indices for the perturbations are 
ns-l = nt - 2(52 = -2 ((5i + $2 + ei) . 



(31) 



(32) 



and the tensor to scalar ratio is r = —8 nt- 

In the SR regime the dynamics of the scalar field can be 

approximated as follows 
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Let us note, from Eqs. (|33l34p . that the Si are not inde- 
pendent of 7 and in the 7^1 limit ^i ^ 7 and thus 
terms proportional to Sl/j ~ Si should be kept in first 
order calculations. Second order terms in the above ex- 
pressions are retained in order to better interpolate the 
regime from large to small 7. On using the above ex- 
pressions and keeping the first order contributions one is 
finally led to [H 



1 = 



27^^ 
1 + 67 



Vcff.crcr* 3Kfr , 
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and 
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7(t; 
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1 + 67 K' 



(36) 



(37) 



where terms proportional to V^^ a*/^* ^^'^ "^^ ^^^^ order 
for 7 <C 1 and second order for 7^1. 



IV. CONSTRAINTS ON DIFFERENT 
POTENTIALS 

On using the above expressions one can investigate the 
observational constraints IG inflation imposes on differ- 
ent potentials for the scalar field. In particular we tested 
the following symmetry breaking potentials which fix the 
Planck mass after infiation: a Landau-Ginzburg (LG) 
potential 



M 



VLG(a) = ^(a^-a, 



a Coleman- Weinberg (CW) type potential 

,4 



Vcw (f 
a cosine potential (CO) 



M 4 / 1 o- 



1 
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Vco (^) = A 



1 + cos I TT 

(^0 



(38) 



(39) 



(40) 



and some generalization of the above potentials of the 
form 



y^(-)-^i-'--iy 



and 



V2 (ct) = A 



n -,2 

^] -l] 



(41) 



(42) 



with n > 1. In particular for our analysis we chose n — 2 
in the case of (PT|) and n = 3/2 and n = 5/2 in the case 

of dUD. 

These potentials share a few interesting features: they 
have a minimum for a = aQ and a relative maximum in 
a = 0, V (cr) > and all of them allow both small and 
large field inflation. We further note that we have in- 
cluded the CO potential even if it agrees to lowest order 
with the LG potential on expanding around cro. Nonethe- 
less it leads to slightly different results because of higher 
order effects. On assuming that the pivot mode exits the 
horizon from 50 to 70 e- folds (A^*) before inflation ends, 
one can compare theoretical predictions for n^ — 1 and r 
with WMAP5+BA0+SN observations [H] for different 
choices of 7. 

In Figures (|2l4p (where the outer border of the lighter 
grey region represents the 68% confidence level and the 
border of the dark grey region the 95% confidence level) 
wc plotted the trajectories of {ns,r) for A^* = 50 (line 
with empty markers) and N^ = 70 (line with filled mark- 
ers) as a function of 7. The markers on the trajecto- 
ries represents particular choices of L^ = log 7 and each 
trajectory splits into a dotted line, for the SF regime, 
and a solid line for the LF regime predictions. Note 
that the transition between the two regimes simply oc- 
curs in the 7 — >■ limit. The dashed line in each figure 



represents the exact consistency condition (1251) which is 
rus = 3r -I- 16 rig — 16. 

Potentials (|38l39p do not constrain 7 in the LF regime 
but they need L^ < —4 for successful inflation in the SF 
case. Indeed from Fig. ([2]) one can see that the markers 
for L^ ~ 3 lie outside the 68% confidence lever region. 
For different choices of the potential wc observe that both 
the SF and the LF regimes impose limits on 7 (apart from 
the potential in (|411 with ji = 2) which is not compatible 
with observations independently of 7). In particular for 
50 < TV* < 70, L^ < -4 is needed in the SF case and 
L^ < —2 in the LF case. 

One can estimate the coordinate of some relevant 
points of the trajectories in Figs. (J2I4|) as functions of 
N^,. On observing than the dynamics of a during in- 
flation is strongly dependent on the value of 7 ([551 [M|) 
then, in the 7 <C 1 limit, one can safely assume that 
|ct* — (To|/(To ^ 1. In order to relate cr* and N^ one 
needs to integrate the SR condition ([55]) by parts: 



1 -H 7 (n(cr) 4- 2) 
7(T (4 — n{a)) 



dcr 



N, 



dN = N, 



(43) 



where n^a) = dlog V^(cr)/dlogCT and we approximate 
(Tend (the valuc of a when inflation ends) with (Tq since 
CTond = 170(1 ± 0(7)). In general, the integral on the l.h.s 
of Eq. (|43)) cannot be solved exactly for any potential 
V{a), and even when it is possible, one further needs to 
invert the result in order to obtain cr* ~ <j{N^) which 
is often a difficult task. For 7^1, however, simpli- 
fications occur and a double series expansion in 7 and 
(ct* — CTo) /ao leads to good predictions. In particular 
one finds a (N^) ~ ctq (l — 2y'jN^), n^ ~ 1 — 2/N^ and 
r ~ 8/A^* for the potentials (|38l39p and for (gH) with 
n = 3/2 and n = 5/2, in agreement with the numerical 
result [2g|. Potential (|^T|) with n = 2 is quite different 
and leads to a (N^) ~ ctq (l - 2y/2^yN^), n,, ~ 1 - 3/iV* 
and r ~ 16/A^* for 7 <C 1. Thus, depending on the shape 
of the potential, for small 7, predictions can be incom- 
patible with observations for a 50 — 70 e-fold inflation. 
For 7 large inflation takes place in the LF regime and one 
can still invert Eq. (|43|) on assuming cr* ^ ctq and series 
expanding 7"^ and (cr* — ctq) /cr* (around zero). One ob- 
tains ris ~ 1 - 2/iV* and r ~ 12/N^ for potential dSH]) 
and n^ ~ 1 - 1.5/iV* and r ~ 4/iV* for potential (|39l) . 
In both cases predictions agree with observations but the 
results are quite different. 

Furthermore in the Figs. pi5|) we have plotted the ob- 
servational constraints on the potentials coming from the 
amplitude of the scalar perturbations 



Pnik*) = (2.445 ±0.096) x 10~ 



(44) 



derived from WMAP5 + BAO + SN |2J. Th e con- 
straints on /i, where /i = A/ctq in the case of (|40l42p 
and /i = A ctq in the case of (pij) , arc plotted for varying 
L^ and assuming A^* = 50 (dotted lines) and A^* = 70 
(solid lines). 



^ 0.2 




FIG. 2: The trajectories of the vector {us, r) on varying 7 and 
potentials S38139\40}^ . For l[38p and fgg|) the markers repre- 
sent L-y = —5, —4, —3 in the SF and L^ = —7, —3, +1 m the 
LF regimes. For (O^ the markers are for L^ = —5, —4, —3 
both in the SF and in the LF regimes. The order of the 
markers is such that dotted and continuous lines join for 
L~f — ^ — cxj. The same order is followed in Fig. Op 



V. COHERENT OSCILLATIONS 

In general, in order to avoid never-ending inflation, it 
is necessary the inflationary potential has a minimum in 
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FIG. 3: Constraints on fi from the amplitude of the scalar 
perturbations \44V for potentials II38139\4^ respectively, where 
fj, = A/(To for \4(Jj) . The dotted line is for A^* = 50 and the 
solid line is for N, = 70 e-folds. The plots on the left refer to 
the SF case and the ones on the right to the LF. 

which case the accelerated epoch is followed by a stage in 
which the infiaton oscillates about it. In realistic models 
of inflation it is usually assumed that during this phase 
the inflaton decays into ordinary matter. This regime of 
coherent oscillations was studied in the context of EG in 
Ref. [1^ through the use of a time-averaging procedure 
which relied on the fact that the frequency of oscillation 
of the scalar field is much larger than the Hubble param- 
eter. 

We have studied the coherent oscillations of the scalar 
field in the context of IG finding the following solutions 

0: 

(j{t) = ao + ^J^sin(c^ot)+Or^j (45) 
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for a potential which can be approximated 

^2 



V{a) 



-{a^aor 
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around the minimum. The frequency of the oscillations 
is given by: 



Wo = 




(48) 
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FIG. 5: Constraints on fi from the amplitude of the scalar 
perturbations \44^ for potentials O^P with n — 2 and \4^ 
with n = 3/2 and n = 5/2 respectively, where /x = Acrg for 
jTTp and ^ = A/ao for g^. T/ie doited line is for N^ = 50 
and the solid line is for N, — 70 e-folds. The plots on the left 
refer to the SF case and the ones on the right to the LF. 
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A. MSA Method 



The method we employ relics on the assumption that, 
during the coherent osciUations phase, two different time 
scales enter the evolution of the field: the frequency of 
the oscillations around the minimum and the damping 
rate of their amplitude, the former being much bigger 
than the latter. This assumption will be verified at the 
end of the calculations. 

On using Eqs. (|3l4p one can write the second order 
equation which governs the homogeneous evolution of the 
scalar field in the absence of matter. For a generic po- 
tential V (cr) one finds 



FIG. 4: The trajectories of the vector (ns,r) on varying 7 for 
potential (OTp with n = 2 and potential U^ with n — 3/2 
and n — 5/2 respectively. In the top figure, for potential pip 
with n = 2, the markers represent L^ — —5, —3 both in the 
SF and in the LF regime. In the middle figure, for potential 
J^g[ ) and n — 3/2, the markers are L-y — —5, —4, —3 in the SF 
regime and L^ = —6, —4, —2 in the LF regime. In the lowest 
figure, for potential [J^ o,rid n = 5/2, the markers are L-y — 
—5, —4, —3 in the SF regime and L^ — —3.6, —2.8, —2.0, —1.2 
in the LF regime. 
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If the potential has a minimum for a = ao then the equa- 
tion governing the dynamics of (5ct = a — ctq can be ex- 



pandcd around it up to second order in Eg j gq obtaining where 
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(5cr /3(l + 67) 
+ ^V 2^ 



; (5^2 + j^ 
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and for small displacements around the minimum of the 
potential one has 



m^ ^ 2 



^ ' -. m^Sa + -5a^ 



dcr 



2' 



©(Ja^ 



(51) 
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where [n] — [m] . The effect of the friction terms in (jl|) en- 
sures that ba I gq <?; 1 for some time t and consequently 
the second order equation is a good approximation for 
t>t. On examining Eq. ([50)1 we observe that second or- 
der contributions are responsible for the "slow" time scale 
evolution while the "fast" dynamics is that of an har- 
monic oscillator with frequency wq- On explicitly taking 
into account the two time scales, we look for a solution 
in the form 

3a^ba{t,T) = (5fTo(t,T)+eJcri(t,r)+C'(e2) (53) 

with T = et. On replacing Sa/aQ -> e Sa/aQ in ([50]) one 
is finally led to the set of partial differential equations for 
6ao and Sai: 



S<7o + CJgfJcro = 

S(7i + uj^Sai 



-Z „ „ ^7^ — ^ — OCTq + Z 



(54) 
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(55) 



The standard MSA method consists in writing the gen- 
eral solution for (l54l) as 



(5cto = ^* [t) e 



Z Wq i 



(56) 



and determining A (r) by requiring the cancellation of 
secular terms in the next to leading order equation (|55p 
[see [25 for details] . From this procedure one obtains the 
following differential equation for A (r) 
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and on setting A{t) = R (r) e' ^(^) one finds 9 (t) = 6*0 
and 



R{r) = 



Ra 
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and Oo = Roy^3fi/j is the inverse of the "slow" time 
scale. Let us note that the ratio r between the two time 
scales depends on the ratio Rq/ctq and on /, where / is 
a function of 7 having values between and \/l/3. On 
setting CTo = Mp/^ one finds 



1 + 67 i^o 
2 Mp 



<1 



i?o 

m7 



<C 



1 + 67 



(60) 



which is the condition for the MSA method to well ap- 
proximate the dynamics. The general solution for (|50p 
turns out to be 



6a 



c^o 



2/r 



1 + rujQt 



cos (wq t + ^0) ■ 



(61) 



and the constants Rq and ^o ^r6 related to the initial 
conditions of the second order equation ([50)1 . 
The MSA method is generally applied to first order, 
nonetheless we found numerically that our results were 
significantly improved once the second order contribu- 
tions on the r.h.s. of Eq. (j55p were taken into account. 
These terms slowly vary in time and give a second or- 
der shift to the centre of the oscillations which has to be 
added to the r.h.s. of (ICTI) 



A^ 



|Ap 



CTO 



(62) 



and plays a crucial role in the calculations of the next 
sections. In particular one has n^G = 6 fiao and new — 
10 fi (To and m^ = 2/i(To for both LG and CW potentials. 



B. Equation of State 

A first application of the above results is the calcu- 
lation of the averaged equation of state of the scalar 
field during the oscillatory regime. It is known (and can 
be shown with the same technique we are employing) 
that the averaged massive oscillations of a scalar field in 
EG are equivalent to a fiuid with null equation of state 
in agreement with |18j . On the other hand, in IG, for 
the potentials we are considering, oscillations around the 
minimum are still massive but the results are somewhat 
different. 

The definition of energy density and pressure for a scalar 
field in IG can be given [11[ further in the EG form: 



p, = 37(T^i?^ Pa = -27fTo [H + 2^1 ■ (63) 

With the above definitions continuity equation is satis- 
fied. At the end of inflation the expressions on the r.h.s. 
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FIG. 6: An example of the average value of the energy den- 
sity of the scalar field, as given in 116 Si) , for fSSH . 7 = 10, 
jj, = 10 calculated exactly (dotted line) and with MSA in S66\) 
(solid line). The plot shows good agreement between exact and 
approximate solutions. 



of (|63p are oscillating functions and one may define the 
averaged quantities (pa) and (Per) where 



1 rt+T/2 
{A{t))=- / A{t')dt' 



T 



t-T/2 



(64) 



and T = 2t: /loq is the period of oscillation. Further one 
may define the averaged equation of state by 



w 






(65) 



On averaging the expressions (p5)) , using the solution for 
Sa + An and keeping contributions up to the second order 
in R{t)/<To <C 1 one finds 



(p,) = (l + 97)(fc) 



(66) 



and 



{Pa)=- 37+v/67(l + 67)sinLJoi {Sa ) (67) 



where 



leading to 



• 2 f v 

{6a )=4woCTo-— -2 

(l + rujQt) 



37 



1 + 97 



(68) 



(69) 



In Fig. ([6]) we plot the average (po-) calculated numeri- 
cally (dotted line) and analytically as in Eq. (|66| for the 
potential ([55)) with 7 = 10 and // = 10. The plot shows 
that multiple scale predictions agree well with the corre- 
sponding exact quantities. Good agreement between the 
analytical predictions of Eq. (|67p and the numerical solu- 
tions is also evident in Fig. ([7]) for the potential ([55)1 with 



0.005 
0.000 



V -0.005 



1 1 "i-^f 



-0.010 



• ■ ■ ■ ■■ ■ ■ ' . .". ■ ■ ~?~-'^ "^•■.~" ^.'^.T"7..-~'..''"'~ 



50 



100 
t 



150 



FIG. 7: In the above plot examples are given of the average 
pressure versus time for the potential ifg^l) with 7 = 10 and 
/i = 10. The dotted line corresponds to both the average value 
of the pressure of the scalar field as defined m fSSf) calculated 
exactly (numerically) and the plot of JfiT] ), since both coincide 
withm the resolution of the figure. The solid line is the aver- 
age of jfiTp and the dashed line is the average of the pressure 
calculated without considering the correction term f62\) . The 
plot shows good agreement between exact and approximate val- 



7 = 10 and /i = 10. In this plot the dotted line represents 
both the numerical (exact) expression and its analytical 
approximation (|67p . since both coincide within the res- 
olution of the figure, the solid and the dashed lines are 
respectively the average of {Pa) with and without adding 
the improvement (p^ to (pTj) . 

In contrast with EG, in IG the averaged equation of 
state defined by (p5|) does not lead to the correct averaged 
expansion, as can be directly verified on using (|46p . The 
averaged Hubble parameter is simply {H{t)) ~ 2/(3t), 
namely that of a matter dominated Universe, indepen- 
dently of the equation of state (|69|) . This difference is a 
consequence of the fact that in EG one can verify that 

{H< 



not hold. 



■^{EGy 



whereas in IG such a relation docs 



VI. REHEATING 

On assuming that the inflaton field is coupled to ordi- 
nary matter, then at the end of inflation such a scalar 
field can decay into lighter particles losing its energy 
density and eventually reheating the Universe. Such a 
process has been traditionally described through a phe- 
nomenological decay width F which enters the equation 
for the scalar field as an additional friction term: 



a + {iH + T)a + Va-&l[H 



0. 



(70) 



On requiring that the newly added friction term does not 
affect the Einstein equations (the Bianchi identities must 
still hold) one observes that the continuity equation for 



10 



radiation should also be modified as: 



PR = -m (pR + Pr) + r^^ 



(71) 



The system of equations (|70l7ip can be solved numer- 
ically but analytical solutions are needed to estimate the 
reheating temperature of the relativistic particle fluid. 
In order to simulate the energy transfer from the infla- 
ton to radiation and estimate the reheating temperature, 
we need the definition of the energy density and pressure 
of the scalar field. In the presence of relativistic matter 
the expressions for the energy density and pressure of the 
scalar field, Eqs. (|63| . must be modified as 



(72) 
(73) 



Per == 3j(tIH^ - PR 



P, = -27^2 [h + ^H^]-pr- Pr- 



with these definitions the continuity equation for p^^, P^j 



is still automatically conserved and the introduction of 
the decay-rate F leads to 



Pa = -'iH {pa 



Pa) - T&' 



(74) 



Let us note that Eqs. ([7T|) and (|74p are formally the 
same as in EG, but they lead to different predictions. 



The reheating temperature can be estimated on 
considering the evolution of pR and {pa) as Trh ~ Pr 
when PR > (pa)- In EG one can naively observe that 
this happens when 3H ^ T and one finds the well-known 



estimate T, 



(EG) 



RH 



yTMi 



The continuity equations for pR and (pa) can be esti- 
mated by the MSA method. Starting from the modified 
equation for a 



J 



1 



1 + 67 



dV{a) ^V{a) 
da a 



2- + - J |- [2 Via) + (1 + 67) a') + 2 pr] + - 
a a \ Zj 1 



67 



0. 



r 



one can expand to second order in 5a ^ a — ao and con- 

• 2 
sider pR ^^ V (a) ^ 6a and T ^ H ^ 6a/ao- These 

assumptions are reasonable for the regime we are consid- 
ering [27|, nonetheless they will be checked numerically. 
On expanding one finds 



Sa + ujQSa + 



n — 4m'^/ao , 
2(1 + 67) 



Sa 



(To 1 + 67 



-Sa-\ 



Sa 
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uJq Sa"^ + Sa" + 



2pfl 
1 + 67 



(75) 



If the "fast" time evolution is still determined by the 
leading order terms Sa^ + Wq Sa^ = then, on setting 
Sao = A (t) y{t) and requiring the cancellation of singu- 
larities, one obtains the differential equation 



2A 



1 + 67 



■A + 3. 



Ipr + {1 + 6j)A2Ef 



37^0' 



^ = (76) 



where Ep = ^y^ 



^y^- On then keeping the second 



order contributions in the continuity Eq. (|71|) and aver- 
aging over "fast" oscillations of y one then finds 



Pr-TAEf+A, 



lpR+il + 6j)A^EF 



37^0 



PR = (77) 



where {Sa ) — A^Ef- On noting that 



lpR + {l + 6-f)A^EF 



■ija^ 



(H) 



(78) 



and 



{pa)=A^EF{l+9j) 



(79) 



where H and pa are expanded to the second order, one 
is finally led to the system of averaged equations 



PR 



<H)pR 



-{Po 



1 + 97 

{pa) = -3{H){pa)--^{pa 



1 + 67 



(H) = 



IpR+'j^ip.) 



370-5 



(80) 
(81) 

(82) 



d(p^) _ 



With 2i^ = (pa). 

On comparing the above equations with those arising 
in EG for a massive inflaton we observe that averaged 
energy transfer is less efficient because of the factors 
(1 + 97)-! and (1 + 67)-! in (ISOlSTI) . 
In Figs. ([8]) we have plotted the comparison of the en- 
ergy transfer obtained numerically on solving (j70p and 
([7T|) with that calculated by numerically solving the av- 
eraged equations (|80ll82p . The continuous line is 



na^ <^' 



and the dashed line is 

^R 



PR + (Pa) ' 



PR 

PR + (P<t) ' 



(83) 



(84) 
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FIG. 8: In these figures the solid line is for i83p and the dotted 
line is for \84^ where both are calculated as numerical solu- 
tions to Eas. i80\l82\) . Circles and squares represents the same 
quantities evaluated on solving {75|7-?[ ) and then averaging the 
definition \7S^ . The figure on top is for F = 2 ■ 10~^Mp, 
7 — 10~^ and the lower figure is for F = 2 ■ 10~^Mp, 
7 = 10. Both situations are described well by the approxi- 
mation scheme used. Time evolution is expressed in Mp^ 
units. 
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FIG. 9: In the above figures the continuous line is pr (in 
Mp* units) as a function of t (in Mp""*^ units) obtained on 
numerically solving Eqs. \8(M8^) . The squares represent the 
same quantity evaluated by numerically solving \7!^71^ and 
the dotted line is pR calculated from the numerical solution of 
fggjgTp . The figure on top is for T = 2 ■ 10~^Mp, 7 = 10~^ 
and the lower figure is for F = 2 ■ 10"^Mp, 7 = 10. The 
figures show that assumptions leading to the expression Ii90[) 
yield good approximate results independently on the choice of 
7- 



The round dots plot ilg- and the square dots plot Clj^ 
where the averages are calculated numerically from 
the exact numerical solutions by using (|64p . In par- 
ticular the plots are obtained on setting 7crQ = Mp^ 
for r = 2 • lO^^Mp, 7 = 10^2 (figure on top) and 
F = 2 • lO^^Mp, 7 = 10 (lower figure). In both regimes 
Eqs. (|80ll82p describe well the exact dynamics and, 
qualitatively, it appears that particle production is 
enhanced for 7 small (when the EG predictions are 
recovered) . 



A. The Reheating Temperature 

We are interested in estimating the energy density of 
radiation when it starts dominating over the scalar field. 
At the onset of this stage the averaged Friedmann equa- 



tion (|82p can be approximated by 



{H) 



l±67/„ 



1+97 



37^2 



(85) 



and the system of equations (|80l81l85p can be solved an- 
alytically. 

On dcf jning p, = ^{pa) and f - (1 + 67)~'f Eqs. 
()80|8ip take the form 



di 

dpa_ 
dt 



37^2 



PR + T Pa 



= -3, 



Sjal 



■ Pa -'^ Pa 



(86) 



(87) 



which is exactly the form one has in EG for a mini- 
mally coupled, massive scalar field. As a consequence of 
this formal analogy one expects pa ~ pr when \/ipa — 
FMp/ (1 -f 67) or equivalently t ^ t'^ . 



Eq. ((87)) can be solved easily and leads to: 

7^0 1 + 97 r2 



(p.) 



3 (1 + 67)^ r^ p,,„ ( r A 
|_^oexp(^ 2(1+67) *J 



- 1 



where Co is an integration constant and we assumed that 
the energy transfer begins at t = 0. On setting 



r =1+1 hd. 1 + 9^ 

° 3Y(p.)o'(l + 67f 



(89) 
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r'7^0 



5/3 
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one can explicitly obtain the dependence of Co on the 
initial condition for the energy density {pa)o- Further 
Eq. ()86p can also be solved and on setting p^i (0) = one 
obtains: 



/c5(T+6^So -Ij f3e5(i+fw)Bo + 5j -60' c^HT+fw) (3 6o + 8) 



20 (1 + 67)' 
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e2(l + 6T)5p _ 1 
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(90) 



where 



So = 1 + r. 
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HPa)o (1+67)^ 
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and 60 = ^0 ^ 1 • If one now considers 60 a-nd sets 7 ctq 
^ 1 



TMr 
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(p,)i/^V3(l + 67)^ 
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(i?) 



<1 



(92) 



since it is generally assumed that F ^ (^0), where {Hq) 
is the average value of H at the beginning of reheating. 
The expression (|90p then simplifies and becomes 



Pfl 



F2 7cr^ 3e2(i+6T) +5_ 



20(1 + 67) e2(i+6T) — 1 



(93) 



and on evaluating it for t ^ {1 + 67) F ^ one finally ob- 
tains: 



Pi,((l + 67)F-i) 



SF^Mp^ 

(1 + 67)' 



(94) 



On performing the same calculations as above in EG one 
obtains 



Pr 



(fJG).p_n^ 3p2^p2 



(95) 



which is the standard result. In Figure pO)) we plot the 
ratio between the reheating temperatures in IG and in 
EG on varying 7. Such a ratio is obtained on evaluating 



T, 



(IG) 



rch 



T, 



(EG) 



rch 



Pr {t*JG) 

(EG) ,, N 

Pr (t*,EG) 



1/4 



(96) 



where p)^ and p]j are the numerical solutions of the 
exact equations and t^ is the time for which pR (<* ) = 
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FIG. 10: The continuous line represents the ratio i96{) 
evaluated numerically on using the exact equations for the 
time of equality t*. The dotted line is the analytical esti- 
mate of the same temperature ratio obtained from \94\y5i) as 



T, 



(IG) ,rj,(EG) 



/rr^ ^ = (1 + 67) 



-1/2 



Pa (t*), and the comparison is made with the analytical 
prediction obtained from (|94l95p . Again we observe that 
the correct behavior is reproduced by our analytical es- 
timates. 



VII. PREHEATING 

The phcnomcnological study of reheating described in 
the previous section was improved by the application of 
the theory of parametric resonance in an expanding uni- 
verse to the fields coupled to the inflaton [13, [2l| . The 
first stage of the decay of the inflaton - dubbed preheat- 
ing - can be rather efficient and leads to a non-thermal 



13 



distribution of the decay products. The subsequent non- 
linear evolution may end up in thernializing the content 
of the universe. 

In this section we study the effect of the parametric 
resonance for the various types of fluctuations present 
and coupled to the inflaton during its coherent oscilla- 
tions. For the analytic treatment we use Eq. (P5|) for the 
evolution of the background, and check our results by nu- 
merical analysis. For the study of fluctuations during the 
coherent oscillations regime the cosmic time is useful and 
one rewrites, through a suitable change of variable, the 
equations of motion of the diverse fields in a Mathieu-like 
form: 



d{ntf 



[A{t) - 2q{t) sin 2nt] y{t) = (97) 



where J7 is a frequency which may depend on the field 
considered. We call the above equation a Mathieu-like 
equation since in the Mathieu equation A and q do not 
depend on time. On taking into account the expansion of 
the universe A and q depend on time and two principal 
effects arise: q decays in time and fluctuations, charac- 
terized by a comoving momentum k, move in the {q,A) 
plane (see for example [25[). An examination of the {q, A) 
plot of the stability/instability regions is useful in order 
to understand if the (time-)evolution of the modes end in 
a stage in which inflaton oscillations cannot be neglected. 
The broad resonance regime (q 3> 1) has been dubbed 
stochastic resonance in an expanding universe [20[. The 
intermediate (q ^ 1) and narrow {q <^ 1) resonance 
regimes are difficult to treat and the exponential growth 
which one obtains with A and q constant in time can 
be completely washed out by the expansion of the uni- 
verse. Let us remember that the small q region of the 
first instability band is: 

2 2 

l-q-^-Oiq')<A<l + q+'L+Oiq'). (98) 

which is situated around A ~ 1. When q is not large the 
effect of oscillations for y arises mainly if the fluctuation 
spends time in the first resonant band. 



Gravitational Waves 



For rujQt ^ 1 this equation can be rewritten as 



+ [Ah + 2g/i sin (ujot)] hk = 



d{uJot/2)'^ 

where A^ and qh arc time dependent functions 
4fc2 / 



Ah{t) 



9 2' ^f^ n J. 



(101) 



(102) 



and we omit terms of order 1/i^. For roj^t ^ 1 the tra- 
jectory in the (q, A) plane is therefore Ah{t) ex k^qfj (t). 
The wavelengths pass through the first resonance band 
but end in the stability region. 



B. Inflaton Fluctuations 

On considering infiaton mode dynamics for a generic 
potential (with a non trivial minimum) Eq. (|15p can be 
recast in the form of a Mathieu like equation by rescaling 
5ak — w a^ Z5(Jk and keeping contributions of the homo- 
geneous mode fluctuations da to first order. The final 
equation for Sak turns out to be 



—rpr- + '"off .a (t) oak 
at'' 







(103) 



where 



Gff,S (t) = 



2 — , .2 



1- 



9/r 



1 + r loq t 



cos UJ^t 



4^/67(1 + 67) . ; 

sm Wot 

97 



-O 



i2 



(104) 



and / and r are defined in (|59l) . Thus, in the rwot ^ 1 
regime, Eq. (|103p reduces to 



d{uJot)^ 



[Aa + 2qa 1 sin (2a;oi) -|- 2qa 2 sin (wot)] Sak = 

(105) 

where Aa- , qa- 1 and qa- 2 are time dependent functions of 
the form 



We begin by considering gravitational waves. The 
equation governing the dynamics of tensor perturbations 
([T9| . on setting hk = a'^^^ahk, can be cast in the form 
of the Mathieu like equation 



dt2 



+ mcs.t (t) hk^O 



(99) 



and the effective mass is 

fc2 



Wefr,t (t) = Wq 



./' 



a-'ujQ 1 + rojot 



sin Wot 



-'i^ 



(100) 



a2wg 



-1, 9^1 = — 7 

Wot 



qa2 = 



277 
2(1-1-67) 

Wot 



(106) 



and we omit terms of order 1/t^. The equation for the 
rescaled IG inflaton fluctuations therefore exhibits two 
oscillating terms: whereas the one multiplied by q^2 is 
similar to the term appearing for gravitational waves, 
qa- 1 does not depend on 7. On neglecting q^ 2 one has 
the following trajectory in the {q, A) plane: 



^.(t) = l- 



I'^/.t-^ 



= 1 



4 (^)" -''«•<-' 
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which shows how fluctuations end in the first resonance 
band asymptotically. The fluctuation 5dk grows linearly 
in time as shown carefully by the MSA method in the 
Appendix. This means that the inflaton fluctuations os- 
cillate with constant amplitude [9], rather than decay if 
oscillations in the homogeneous inflaton are averaged in 
time. Since q^ i does not depend on 7 our analysis applies 
to the coherent oscillation regime of chaotic quadratic 
inflation in Einstein gravity as well, in agreement with 
Ref. [24I. Therefore during the coherent oscillations of 
a massive inflaton, gauge-invariant inflaton fluctuations 
oscillate with constant amplitude on large scales [231 and 
also on small scales. We have numerically checked that 
the inclusion of the oscillation term multiplied by 9^2 
does not modify the qualitative behaviour just described. 
Indeed the amplitude of this term is never important, 
go-2 ~ 0{'^/{ujQt)) for 7 < 1 and 9^2 ^ 3/(2a;ot) for 
7^1, and the frequency of oscillations is half of the 
dominant one. 



C. Scalar Fields non-minimally coupled to gravity 

We finally take into consideration the evolution of a 
scalar test field x non-minimally coupled to gravity and 
interacting with the inflaton. Let 






Cu 



be the action for such a field and 



^in 



e 



Rx' 



9 2 2 
— cr y 
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the interaction with ^ > 0. On considering only the inter- 
action of the modes of the field x with the homogeneous 
part of the inflaton one finally obtains 



-1^ + TOcff.x (*) Xk = 



(110) 



where 



mcs.x (t) = ^0 



P 



I + rujQt 



sin LjQt 



+2g 



2I + 67 fr 
'I 1 + rujQt 

2 I + 67 



sincjQ^ 
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(111) 



and Xk 



1,3/2 



Xk ■ From (|111[) we see that the scalar field 



does not end up in the instability region unless 

^l I .,2 1 + 67 ,, 1 
^2 ^-^ 2n ~ 2' 



(112) 



which shows how a mass m'^ generically suppresses the 
resonance. Note also that a broad resonance regime is 
possible for tti^ = g = for large ^. 



VIII. CONCLUSIONS 

We have investigated in detail inflation in the IG 
framework where inflation is driven by the same scalar 
field as that associated with the observed value of New- 
ton's constant. Inflation in IG can be successful in the 
SR regime, and leads to a nearly flat spectrum of scalar 
perturbations and a small tensor to scalar ratio compat- 
ibly with observations. The SR conditions are not sim- 
ply associated with the shape of the inflaton potential 
since the dynamics of the scalar field strongly depends 
on its coupling 7 to gravity. We have made accurate 
comparisons with observations for different "symmetry 
breaking" potentials and showed how their shape is rele- 
vant for successful predictions of the model. In particular 
we have seen that the LG and CW type potentials are 
compatible with the data for large field inflation and for 
suitable values of the parameters for small fleld inflation. 
For potentials involving higher powers of the scalar field 
(with respect to the previously mentioned cases) agree- 
ment with the data is more problematic and if it occurs 
one obtains constraints on 7 for both the large and small 
field cases. 

We have also studied, always in IG, the coherent os- 
cillatory regime and obtained an accurate analytical so- 
lution for the dynamics of infiaton and of the Hubble 
parameter by using a MSA. We employed the solution to 
first investigate the behavior of the average evolution of 
the Universe. We found that, independently of the av- 
erage equation of state of the scalar field, the Universe 
expands as {H{t)) oc 2/(3t). 

Further perturbative and resonant reheating have also 
been studied. We found, in the former case, that the de- 
cay of the inflaton into ordinary matter is less efficient 
for 7 large and that the reheating temperature is com- 
parable with the EG value for 7 small. We also exam- 
ined the parametric ampliflcation of the scalar and tensor 
perturbations and of a generic scalar field non-minimally 
coupled to gravity and coupled to the inflaton. We note 
that parametric resonance occurs for scalar perturbations 
with H <ti k/a ^ wq and such perturbations maintain a 
constant amplitude instead of decaying as the Universe 
expands. Finally parametric resonance does not occur 
for gravitational waves whereas amplification may occur 
for a generic scalar field depending on the parameters of 
its Lagrangian. 



Appendix A: Preheating 

In an expanding spacetime, parametric resonance in 
the first instability band occurs when the Mathieu-like 
equation governing modes dynamics can be cast in the 
form 



y+[l + eq{t)sm2t]y = 



(Al) 



where e is a small dimensionless parameter. The theory 
of parametric resonance is well known when q (t) = const 
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but in an expanding spacetime, where a time dependence 
generally appears, the dynamics leads to quite different 
results. MSA gives the correct approximate behavior for 

(ED. 

Since eq (t) <^ 1 one can look for an approximate solution 
in the form 

y^yo{t,T)+eyi{t,T) (A2) 

where t = et. On comparing orders one has 

d'yo 






yo = 0=^ya = A{t) e** + c.c. 



2/1 



dtdi 



qit)' 



2i 



(A3) 
yo = (A4) 



and the requirement that secular contributions cancel in 
^4p leads to the differential equation 



d4 

4-_,(t)A*.0. 



(A5) 



On separating A into real and imaginary parts {A — B - 
i C) growing and decaying modes also separate 



4f + ,(0C 



(A6) 
(A7) 



and, if q (t) > 0, the imaginary part decays. 
In general, the expansion of the universe during reheat- 
ing leads to q (t) ^ p/t" and in this case the growing 
dynamics becomes 



B (t) = Bo exp 



to 



At"- 



(A8) 



When n > 1 the amplitude of oscillations increases 
but still exhibits a constant asymptotic behaviour B ^ 



So exp 



.t' 



for n < 1 the amplitude increases 



4(1- 



,t' 



In contrast. 



4(n-l)''0 

exponentially as B ^ Bq exp 

when n = 1 (which is the value we find for scalar pertur- 
bations in the instability region), one obtains the follow- 
ing power-law dependence 



B = Bo 



to 



p/4 



which, for p = 4 corresponds to 5a 
Section VIIB. 



(A9) 
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